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We propose a method to reconstruct the vibrational quantum state of molecules excited by a 
general excitation laser pulse. Unlike existing methods, we do not require the molecules before 
excitation to be in a pure state, allowing us to treat the important case of initially thermally excited 
molecules. Even if only a single initial level is appreciably populated, initial levels with small 
populations can still give major contributions to the unknown vibrational state, making it essential 
to take them into account. In addition to the excitation pulse, the method uses two incident, 
short laser pulses in a non-co- linear geometry to create four- wave mixing in the molecules. The 
measurements used in the reconstruction are spectra of the outgoing four-wave mixing pulse at 
different time delays of the excitation laser pulse. An important point is that the method does not 
require detailed knowledge of molecular transition moments between excited states nor of any of the 
incoming laser pulses, but circumvents this requirement by using one or more calibration laser pulses 
in a separate experiment either before or after the main data are recorded. The only requirements 
for the calibration laser pulses are that the constant parts of their spectrums should together cover 
the spectral range of the excitation laser pulse, and the constant part of each should have sufficient 
spectral overlap with one other calibration pulse to populate two of the same levels. Finally, we 
discuss the extension of the reconstruction method in this paper to more general situations, hereby 
presenting the new idea of quantum state reconstruction through perturbations with calibration. 

PACS numbers: 03.65. Wj, 42.65.-k, 42.65.Dr 



I. INTRODUCTION 

Everything worth knowing about a molecule at any 
given time is its physical state. Hence, whether one's in- 
terest lies in interaction of radiation with the molecule, 
chemical reactivity or intra-molecular processes, the state 
contains this information. Evidently, it is of great inter- 
est to be able to determine the state of a molecule. In 
this theoretical paper, we show how one can determine a 
molecular state in a four- wave mixing experiment. The 
method presented is particularly useful for determining 
the molecular states prepared by a single, optically tai- 
lored femtosecond laser pulse. Production of this type 
of laser pulses has been demonstrated and they 

have been used to produce specific vibrational states in 
molecules 0- 

Finding the state within a quantum mechanical de- 
scription is complicated, as there is no single observable 
giving the quantum state. This is well-known for pure 
states, i.e. states that can be characterized by a single 
wave function, where there is no observable directly re- 
vealing this complex wave function. The same absence is 
true for the more general case of mixed states, which are 
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states that can be characterized by a density matrix, but 
not necessarily by a single wave function. In this paper, 
we shall be treating such mixed states, paying special 
attention to thermal states; a type of mixed state that 
describes a quantum system at a finite temperature. In 
particular, we will concentrate on finding the vibrational 
state of an initially thermal molecule, which has been 
excited to an unknown vibrational state by a short laser 
pulse. 

Although the state is not directly observable, one can 
determine the state from measurements. This is the field 
of study known as quantum state reconstruction. On 
each member of an ensemble, one or more commuting 
measurements are performed, whereafter this member is 
discarded. This ensemble must be large enough that the 
measured values are close to their true expectation val- 
ues, which obviously requires numerous measurements 
of each quantity. Fortunately, in the method presented 
below there will be a vast number of molecules in each 
experimental run, ensuring excellent statistics. 

In the broader picture, there has recently been grow- 
ing interest in reconstructing quantum states for various 
systems, for a review see [6[. The general problem of 
quantum state reconstruction is experimentally and the- 
oretically challenging, and methods tend to apply to only 
a very particular setting. The most well-known is the 
tomographic reconstruction of harmonic oscillator states 
from measurements of spatial distributions Q , [H • Other 
systems include particles in traps (neutral atoms Q and 
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ions [HI), general one-dimensional systems [Tl|, dissoci- 
ating diatomic molecules [12] and the angular state of a 
dynamically aligned molecule [HI . 

Earlier works directly related to this paper have 
proposed reconstructing vibrational states for diatomic 
molecules using heterodyne detection of fluorescence [13] , 
for general molecules using time- and frequency- resolved 
fluorescence [Hj], [Tf3], and using time- and frequency- 
integrated fluorescence and a known reference state for 
quantum state holography p7| . Common for these pro- 
posals is the requirement that the initial state be a pure 
state. Realistically, fulfilling this requirement in turn de- 
mands a sufficiently low temperature of the molecules 
before the creation of the vibrational excitation. This 
initial temperature must be so low that only the lowest 
vibrational level is populated. For many molecules this 
is unpractical, e.g. gaseous iodine, where practical tem- 
peratures populates more than one vibrational state, lest 
the iodine is deposited [Hj]. It is important to note 
that even though the initial thermal state's ground vibra- 
tional state population is usually much larger than that 
of all other vibrational states, it does not guarantee that 
this gives the dominating contribution to the unknown 
excited state. For instance, in the experiment on molec- 
ular iodine described in [l9| , the ground and first excited 
vibrational level in the initial thermal state give approxi- 
mately equal contributions to the unknown excited state, 
even though the ratio of populations is about 2.5 : 1. This 
is mainly due to the different Franck-Condon factors in- 
volved in the formation of the unknown excited state by 
the first laser pulse. Hence, taking several initially ther- 
mally populated levels into account can be essential to 
understanding the excited state formed. 

The paper is arranged as follows: In section [II] we iden- 
tify the type of quantum state being reconstructed and 
give a brief description of what is required of the mea- 
surements. In section ITO1 we outline an experimental im- 
plementation of the required measurements and give an 
overview of the energy- and time-regimes involved. In 
section IIVI we show by calculation how to reconstruct 
the unknown vibrational quantum state. In section [V] 
we discuss the generality of the reconstruction procedure 
and conclude the paper. 



II. STATES AND MEASUREMENTS 

Having established the importance of the physical state 
and its reconstruction, we will now specify what we mean 
by the physical state and what the measurements for re- 
construction must fulfill. 

It is essential to distinguish between the state of a sys- 
tem and its dynamics. The state contains the answer to 
all questions we may ask the system at a certain point 
in time, i.e. the result of all thinkable one-time measure- 
ments at this time. In sharp contrast to this state is the 
concept of dynamics, which is prescriptions for finding 
the state at a later time, from knowing it at an earlier 



time. 

In the subject of state reconstruction, and in this pa- 
per in particular, the goal is to determine the state of a 
system. The type of state we will find in this paper is 
an unknown vibrational state of a sample of molecules, 
coherently excited from an initial thermal state. Specifi- 
cally, we will restrict the reconstruction of the unknown 
vibrational state to the electronic state(s) it populates, 
denoted by the electronic index a, and thus not find cor- 
relations with the initial thermal population. This is ex- 
actly what is usually meant by the concept of an excited 
vibrational state. For transparency, we restrict ourselves 
to the most important case where the initial thermal state 
is limited to the electronic ground state, with possibly 
several vibrational levels populated. The excitation that 
forms the unknown vibrational state need not be pertur- 
bative in the sense that it may transfer a large fraction 
of the initial thermally populated level to the electronic 
state a. Nevertheless, it must leave a non-vanishing pop- 
ulation in the initial state and not transfer population 
back to other of the initially populated levels. 

Since the initial state is a mixed state, we cannot de- 
scribe it by a single wave function, and because coherent 
excitations preserve the mixed character, neither can we 
describe the excited vibrational state this way [2(j. A 
sufficing description is the density matrix, for which the 
initial thermal state would be diagonal in the energy ba- 
sis, but we shall take a more transparent course. 

As the states originating in different energy levels in 
the initial thermal state are incoherent with respect to 
each other, and the vibrational state we seek is coher- 
ently excited from these, we can describe the system as 
a incoherent sum of vector states (corresponding to wave 
functions). By an incoherent sum, we mean forming a 
sum of the states, each with a phase factor exp(i6k) la- 
beled after the original energy level k in the initial ther- 
mal state from which the coherent excitation took place. 
These phase factors will be averaged over when taking in- 
ner products, making incoherent cross-terms vanish. In 
this paper, each incoherent term will be labeled by a left 
subscript, such as To illustrate this, consider the 
example of a thermal state of a molecule with the only 
energy levels populated being the vibrational levels k = 
and k = 1 in the electronic ground state. We then imag- 
ine exciting a vibrational state in another electronic level 
a, the vibrational eigenstates \<p a ,i) herein being labeled 
by index I. The excited state will then be: 

IV>a) = ^2 \k1pa) 
k 

= |fe=o^a) e i6 ° + |fe=l^a) e** 1 

= e 4eo ^ o A|^)+e l9l ^iA|0a,i), 
i i 

where the are expansion coefficients. What we mean 
by "finding the unknown vibrational state" is thus deter- 
mining all the complex expansion coefficients kPi- With 
the incoherent averaging, inner products between states 



originating from different levels in the initial thermal 
state vanish, e.g. 

( k =^a\A\ k=1 ^ a ) cx fie- 9 " /^e'^O, 
Jo Jo 

where A can be any operator. Thus, we will not see any 
cross-terms between states originating in different energy 
levels of the initial thermal state, whence we will conve- 
niently treat these states for different k separately. How- 
ever, one should notice that it is still possible to see in- 
terference in an intensity signal between the light formed 
from different incoherent levels. This can be easily seen 
by considering the two electric fields formed from the 
transitions from level a to the electronic ground state 0: 

k=oE(t) cx d 2 1 dt 2 (k=otpo\d\k=oipa) 
k=\E{t) cx d 2 /dt 2 ( k= iipo\d\k=iipa}, 

where d is the scalar product of the dipole moment oper- 
ator with the electric field's polarization vector. Since 
the incoherent 6>fc-factors cancel in each electric field, 
we will observe interference in the intensity signal / oc 
[ k =oE(t) + k=iE(t)f . 

Having accounted for what we mean by the state, we 
now turn our attention to measurements. To accomplish 
the state reconstruction, we will require knowledge of the 
results of enough measurements to uniquely determine 
the state [2l[ . Such a set of measurements is known as a 
quorum. 

Finding a quorum can be straightforward or tedious, 
depending on both the system and the level of descrip- 
tion. In a classical system of particles, one quorum is 
given by measurements of the coordinates and momenta 
of all particles at a fixed point of time. In contrast, the 
quantum mechanical description used in this paper makes 
finding a practically realizable quorum much more chal- 
lenging. We will use the prevalent approach which is to 
perform the same type of measurement at different points 
of time, thereby letting the dynamics reveal the state, see 
e -g- 0i [HI- Using this approach, one is almost always 
forced to assume full knowledge of the dynamical laws of 
the physical system. However, we shall largely circum- 
vent such assumptions by instead performing a calibra- 
tion of the measurement apparatus. Specifically, we will 
show that spectra of a four- wave mixing laser pulse, cre- 
ated at different times after the formation of an excited 
vibrational state, encompass a quorum for determining 
the set {k0i}. 

III. EXPERIMENTAL IMPLEMENTATION 
OUTLINE 

In this section we outline a pump-probe experiment 
where one can record the data necessary for the vibra- 
tional state reconstruction. This is a four-wave mixing 
experiment where one sends in three laser pulses 1-3 from 
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different directions on a molecular sample, and records 
the spectrum of the outgoing four-wave mixing (FW) 
pulse , [3 • Laser pulse 1 reaches the molecular sam- 
ple first and can be regarded as a pump pulse, since it 
has the effect of creating the unknown vibrational quan- 
tum state in the molecules. Pulses 2 and 3 arrive at 
the molecules long after pulse 1 is over, and can thus 
be regarded as probe pulses . The time r when pulse 1 
creates the unknown state is varied from one experimen- 
tal run to the next, retaining the timing of pulses 2 and 
3. One then measures the spectrum of the resulting FW 
pulse. The time sequence of the four laser pulses is out- 
lined in Fig. [TJ One should notice that even though the 
pulses 1 — 3 should individually be the same from one 
run to the next, there is no requirement of relative phase 
stabilization. 




Figure 1: Time sequence of the laser pulses. First, pulse 1 
excites a vibrational state in a sample of thermal molecules at 
time r. When pulse 1 is over, the molecules are probed using 
the pulses 2 and 3. While the pulses 1-3 are always the same, 
the time r is varied from one experimental run to the next. 
The interaction of the three pulses 1-3 with the molecules 
gives rise to a four- wave mixing (FW) pulse. The spectrum 
of this FW pulse is measured together with the current value 
of t. Of further notice is that the emission of the coherent 
FW pulse can continue long after the incident pulses 1 — 3 
have died out, since the emission is due to an excitation in the 
molecules. For transparency, we will in section llVl assume that 
pulse 2 and 3 have no temporal overlap, as shown. However, 
this is not a requirement of our method, and we will relax 
this assumption in section [V] allowing for overlapping probe 
pulses as used in [jg| |. 

To keep the theoretical treatment transparent, we will 
now introduce several simplifying assumptions. We em- 
phasize that these are not formal requirements, and are 
only included for clarity. Indeed, we will do away with 
them in section IVl 

The laser pulses involved will all cause transitions in 
the molecules, dominated by dipole transitions. In Fig. [2] 
we show an overview of such a series of transitions, each 
being accompanied by a change of the electronic state. 
We have sketched the four electronic states and a-c, and 
within each of these several vibrational states. Further- 
more, the figure shows that the transitions due to pulse 2 
always occurs before transitions due to pulse 3. Realisti- 
cally, this could be the case if pulse 2 precedes pulse 3 so 
that they have no temporal overlap, as shown in Fig. [T] 
We also require that only a single photon from pulses 2 
and 3 are involved in each transition. Finally, we neglect 
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the rotational degrees of freedom. This constitutes the 
most transparent situation, and we shall initially calcu- 
late the spectrum of the FW pulse in section IIVI using 
these assumptions. In section fVl we will then discuss how 
to relax these non-essential assumptions. 

One can experimentally justify the implicit assump- 
tion made in Fig. [2] that only a single photon from each 
of the incoming laser pulses is involved in each tran- 
sition. Using the non-co-linear, so-called folded BOX- 
configuration shown in Fig. [3j one sends in pulses 1-3 
from three different directions. Each angle between the 
pulse propagation direction is small, typically a few de- 
grees. One can select out the FW signal caused by pre- 
cisely one-photon interactions from each pulse, simply 
by using pinholes to select out the appropriate phase- 
matched direction for the FW pulse. 

An important point, which will play a key role in the 
calculations below, is that the directional selection will 
typically only select the number of photons involved in 
each transition, and thereby the electronic states, but 
not select which vibrational levels are involved in each 
transition. The reason is that the electronic states are 
separated by a comparatively large energy, giving rise to 
a realistically observable difference in angle for the FW 
pulse. In contrast, the vibrational levels within a certain 
electronic level are energetically so close that the focus- 
ing angles of the incoming pulses are much greater than 
the change in angle needed to satisfy phase matching to 
one vibrational state or the other. The ensuing interfer- 
ence in the FW intensity is precisely the phase-sensitive 
quantities that we shall use in the reconstruction. 

A hnal important reason to use the non-co-linear ge- 
ometry is that it allows one to observe the FW pulse 
without the presence of the very powerful background 
from pulses 1-3. 



the time shown in Fig. [TJ the time when the unknown ex- 
cited state is formed. This state ^^(t)) does not need 



to be formed perturbatively, as long as it leaves a non- 



k^pQ 1 ^ -state, and it does not 



vanishing population in the 

transfer population back to other of the thermally pop- 
ulated levels (i.e. other fc-indices). We can then expand 

the unknown state k' t Pc^\ T )^ on t ne vibrational eigen- 
states \(f> a ,i) with energy fuvf, belonging to the electronic 
level a: 



^kPllta.de^e-^, (2) 



where the kPi are expansion coefficients and the transi- 
tion frequencies between the vibrational states v 
and v' in their respective electronic states a and a 1 are 
given by: 



Unlike pulse 1, the probing pulses 2 and 3 must be per- 
turbative. Using the electric field E-iit) of pulse 2, we 
find: 



| Idle h E,U-2) 



x d 



(3) 



where H a is the vibrational Hamiltonian in the electronic 
state a, q = — e is the electron's charge, and d is the 
dipole moment operator. The corresponding expression 

for 



ktyc (t) ) can be found similarly: 



IV. THEORY 

In this section, we shall theoretically treat the situation 
outlined in the previous section, aiming at reconstructing 
the quantum state formed by pulse 1. We will accomplish 
this through calculating the spectrum of the FW pulse. 

We recall from section |TT] that we initially have a ther- 
mal state, where each energy level k may be treated in- 
dependently on the quantum level. Concentrating on a 
fixed k, the total state formed by the three laser pulses 
can be written as: 

*V>(r)) ( = \k4 0) ) t + \k^Kr)) t + 

k4 2) (r)) t + \ k ^Hr)) t , (1) 

where the superscript denotes which pulse has formed the 
state from the previous, the right subscript denotes the 
electronic state, and t is the time. In each of the pertur- 
bation orders, we shall for clarity retain only terms lead- 
ing to coherent FW emission. The parameter r signifies 



K 2 



dU 



dt2e -iH c {t-t 3 )/h E ^ 



OO J — CO 



x de- l ""^- t2 ' )/h E 2 (t 2 ) 



x d 



1-2 



(4) 



We introduce the Fourier transforms of the electric fields: 

/oo 
duj r £ r {cj r )e- llAJrU , r = 2,3. (5) 
-oo 

Furthermore, we use identities resolved on the vibrational 
states v of the electronic state a, including both bound 
and continuum states: 



la = \<t>a,v) {(t>a,y\ 



(6) 



Finally, we abbreviate q/H times the dipole transition 
moments between the vibrational levels v and v' in their 
respective electronic states a and a': 



y) 



(7) 
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Figure 2: Schematic energy diagram showing the levels of the molecular system and the four laser pulses with time-dependent 
electric fields Ej(t), j — 1,2,3, FW. Assuming the molecules are initially thermally excited, we can consider each initially 
populated level k in turn, and perform a incoherent sum in the end. The system starts out in the electronic state 0, wherefrom 
laser pulse 1 creates a vibrational state in the electronic state a. After a time r, which is varied from one experiment to the 
next, the state created with pulse 1 is probed through a four-wave mixing process using pulses 2 and 3 to transfer the system 
through electronic states b and c. Under coherent emission of a FW photon, the system finally ends up in the k where it 
originated. Some of these electronic states may be identical, e.g. = 6 and a — c in the experiment described in The 
delay time r is chosen small enough for pulse 1 to be over when pulse 2 and 3 are applied. 



We use Eqs. ©-(0) with Eq. Q and perform the time 
integrals. To accomplish this, we make use of the usual 
trick of letting £ ' a , and (D^'° have a small negative imag- 
inary part, whereby the integrals converge. 

poo poo 

1, 1 _ J -co J-oo 



k,l, 



^2(^2)^3(^3) e 



-1(^2+0)3+1!" )t 



(. ~ a.b \ ( ~ c.a 



LU 2 - W3 



One may note in passing that if the dominating effect of 
pulse 2 is to stimulate emission as shown in Fig. [2j the 
negative frequency components of £2(^2) will give the 
dominating contribution to the ^-integral. 

The electric field of the FW pulse is proportional to 
the second time derivative of the dipole moment 



E FW (t,r) = 22kE FW {t,T) 

k 

d 2 



This yields for the contribution caused by the fc'th ini- 
tially populated level: 



^m(^)] = ^^mwn 



-,<°)\ 



(8) 



Notice, that we only have to use the same fc'th vibrational 
state from the initial thermal state, since inner products 
with all other terms vanish (see section [TTJ) . Using this 
result, we can find the frequency dependence of the FW 
pulse's electric field: 

k£Fw(u,r) 



1 

27 



dt e iuJt 



4{t,T) 



J2kC{(u;)eOk[3I 



(9) 



where the mapping coefficients kCi{ui) (and their com- 
plex conjugate) depend upon the lifetime of the energy 
levels and the detailed electric fields of the pulses 2 and 3: 



I ,m ,n M n ,k 



00 £2(^2) ££ f^fc '? - &2 - Uj) 

dcvo ^— '- 



U) 2 



~ a,b 
■ u t'm 



(10) 



This expression clearly shows the transitions a — > b — > 
c — > through the ©-factors and the resonance frequen- 
cies in the w-dependent denominators. The apparent di- 
vergencies are not real, but a usual artifact of applying 
perturbation theory. 

The most important point for our purposes are the 
linearity of Eq. (0 in our sought coefficients kfl?, an d 
the exponential r-dependence. 
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Figure 3: Geometric outline of the laser pulses' interactions 
with the molecular sample. Dashed lines signify propaga- 
tion below the plane of the paper, while fully drawn lines 
signify propagation above. It is central to the calculations 
in this paper that one uses the shown non-co-linear config- 
uration, in which the pulses 1-3 have different propagation 
vectors (BOXCARS-configuration for pulses I and 3 identi- 
cal). As seen by the molecules, this will be similar to the 
pulses 1-3 approaching from three of the four corners in a 
rectangle. By using a pinhole, one may separate out the FW 
pulse, hereby ensuring that the measured interaction in the 
molecular sample has involved precisely one photon from each 
of the pulses 1-3. It should be noticed that, in contrast to 
usual four-wave mixing situations, the BOX condition gives 
only state-selectivity of the electronic state, but not of the 
vibrational states herein. This practical inability to tell apart 
vibrational levels by the angle of the FW photon arises due to 
the energetic closeness of the vibrational levels compared to 
the energetic distance between the different electronic states. 
Thus, transition to a certain electronic state gives rise to a 
easily detectable angular difference, whereas the transition to 
different vibrational states herein does not; see also Fig. [2] 



The measurements at our disposal are spectral inten- 
sities for each delay r, i.e. a two-dimensional spectro- 
gram. Even though we measure intensities and not elec- 
tric fields, the signal is found from squaring the electric 
field S(u,t) oc \£ fw (lu,t)\ 2 [27]. Therefore, dividing the 
signal by lo a , we can find the following quantities: 



the chosen frequency f2: 
1 



F(Q,uj) = 



2tt 



N— I (It 

2tt 



d(u,T) 



-■zOt 



N^'kCiMkftvCpWvfi, (12) 



with the primed sum running over indices (k, k' , I, V) so 
that VL = ui^'p — Even though the r-integral in 

Eq. (fT2|) runs over all values of r, one cannot achieve this 
in real experiments. Both the inevitable truncation of the 
r-interval and decay processes will cause a broadening of 
the peaks in F(Cl, ui) as a function of f2, in turn leading to 
an increased number of frequency recurrences. Although 
a concern in principle, it will be seen that it is possible in 
typical experiments to use (a discrete version of) Eq. Q12j) 
to find sums with just a single or a few terms. We will 
aim to extract single terms from these sums. Obviously, 
choosing Q = gives many terms, among them the ones 
with 1 = 1', making it difficult to extract any single one. 
Fortunately, it will turn out that we will not need these, 
so we shall henceforth look for terms with I ^ V . 

To select out a single or a few terms from the function 
F(Q, ui) in Eq. (|12p we must make appropriate choices of 
the variables O and ui. In the following subsection IIVA1 
we shall present a method to do this for the k that can be 
selected by proper choice of ui: case (I). This is usually 
the case for the lowest populated (ground) state k = 
(see e.g. [l||) and for any k that are energetically sep- 
arated from their neighboring vibrational levels by more 
than the combined width of the laser pulses. Further be- 
low, in subsection IIV CI we will present how to deal with 
the remaining k through making a weak assumption on 
the spectrum - case (II) . Finally, in subsection IIVDI wc 
will discuss how to do the reconstruction if the two above 
methods fail. 



= NY,kCi(uj) k ,Cf,(Lu) k f3i k ,/3t, 
i,i' 

k,k' 

xe i ( ffl r.i?- fi> 2:°*') r j (ii) 

where the overall positive factor N is introduced for two 
reasons: First, we do not know the absolute popula- 
tions created by the pulses, and second, we wish to avoid 
the difficulties connected with measuring absolute signal 
strengths. To extract the kfiis, a first step is to iso- 
late single terms kCi{u>) kfli k' Cf, (uj) & (3*/ from the sum 
in Eq. We can use Fourier transformation with re- 

spect to t to extract a single or a few terms of this type. 
The number of terms found from such a procedure equals 
the number of recurrences in the transition frequencies at 



A. Emission due to single fc's: Case (I) 

We turn our focus to the group of k states that belong 
to case (I). Choosing the frequency fl = uifft, we see that 
the sum Eq. (fl2| contains terms with all k — k' . Even 
though this is a complication, we can circumvent it in 
this case by using the freedom in choosing ui to select k 
by making all k Ci(ui) vanish, except for a single k. This 
situation is illustrated in Fig. 0] 

In the usual experimental situation where the signal 
S(uj,t) is simultaneously recorded for all frequencies cu 
in a wide range, we also know F(Cl, ui) for this range of 
ui. To find out if only one k contributes at a fixed ui, 
one simply checks whether F(ujf)f — w fe ' fc; ,cj) = for all 
k' ^k. 

To avoid testing too many values of u> by the above 
procedure, one can readily determine the approximate 
u-range of interest. From the easy measurements of the 
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central frequency and widths of pulses 1-3, one can locate 
a limited w-interval where one expects the FW signal 
from each k. In such considerations, one can also ben- 
efit from knowing approximate potential surfaces to get 
an idea of which transitions have appreciable transition 
matrix elements. 

Thus, by wisely selecting certain values of w, we can 
limit the primed sum Eq. (|12p to run over only I and V . 
If there are degeneracies in the frequencies u)f\?, it can 
still be possible to find the elements indirectly. Since the 
details depends on the problem at hand, we illustrate this 
with an example. 

Frequency uj 



ui 



k = 1 dominates 



Both k = and k = 1: 
Interference 



k = dominates 



can be difficult to measure NCi fti C\, /3f, with the integers 
I and I' differing by more than a few. This is because 
terms with I and I' significantly different correspond to 
quantum beats between distant levels in the unknown ex- 
cited state. Fortunately, these are not necessary because 
of our freedom of choosing normalization and phase of 



each 



as will be clear below. 

We have now reached the limit of what we can ex- 
tract from the (oj, r)-spectrograms. To find the unknown 
quantum state through the Pi 's, we need to know the 
Cj(o;)'s. These are in turn both difficult to calculate, 
one reason being the regularization of the apparent di- 
vergencies in Eq. (|10p . and more importantly require 
detailed knowledge of the electric fields of the pulses 2 
and 3, of the dipole matrix elements 2)"'", and of the 

transition frequencies to"'", ■ One may rightly argue that 
with the knowledge of the dipole transition moments, it 
would have been easier to measure the field of the single 
pulse 1 and simply calculate the unknown quantum state 
by propagation using the dynamical laws. Therefore, we 
shall take a completely different approach which circum- 
vents both regularizations, detailed knowledge of electric 
fields, and knowledge of transition matrix elements and 
-energies between excited states. 



Figure 4: Schematic illustration of a spectrogram Eq. 
in case (I): There are regions of uj where only a single k 
contributes. The vertical lines signify maxima in the signal 
S(ui,t)/u> , with contributions from different initial levels ex- 
plicitly shown. To extract the numbers kPi from the sum 
Eq. 1)120 . one can choose uj to select k. In the example shown 
in this figure, one can extract information about of5i from 
using the data with uj — u>o- Similarly, fixing the frequency 
at u) = u)i allows one to obtain the coefficients i/3;. 



Say, uj 



~ a, a 



and we have Fourier transformed 



to find the following sum, suppressing the fc-index and 
w-variable for clarity: 

F(Q) = NQ^Cp^+NQ^Cp^. (13) 

Provided that we can find each of the following coeffi- 
cients and that none of them are zero, we can form the 
quotient: 

NC - NC^C^,, 

allowing us to extract each of the terms in Eq. (fT3"|) . Thus, 
we can use intermediate levels to obtain each term in 
sums, otherwise hidden because of the degeneracy in the 
transition energies. 

In principle, we could similarly find all such products 
NCiPiC*,/3*,, including 1 = 1'. Though having all these 
products permits a nice statistical treatment through sin- 
gular value decomposition of the outer product matrix, it 



B. Calibration pulse 

In this subsection, we introduce one or more calibra- 
tion pulses aimed at creating excited quantum states in 
the electronic level a with known calibration expansion 
coefficients k f3f al . This is a different approach than het- 
erodyne detection or other reference pulse techniques, 
since the calibration experiment is a separate experi- 
ment from the one where the reconstruction data are 
recorded. Hence, by performing exactly the same exper- 
iments as above, except with the calibration pulse sub- 
stituted for pulse 1, we will find products of the form 
Nc kCi{u>) k'C*,(uj). These will in turn be used in the 
reconstruction experiment with pulse 1 to find the coef- 
ficients fc/3/. Let us consider this in detail. 

Contrary to pulse 1, we will require that the calibra- 
tion pulse(s) be perturbative, because this gives simple 
expressions for the state created in the calibration. We 
find this state similarly to Eq. letting the upper limit 
in the integration tend to infinity since the electric field 
E ca i(t' + t) vanishes for times long before the pulses 2 
and 3 arc turned on: 



x d e -**o(*'-o)/» 



fc-00 



(0) 



x e 



-iujft2 



(14) 
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Comparing Eq. (fT4f to Eq. ([2]) , we see that the expansion 
coefficients of the calibration state are 



o,0 c f~a.0\ 
I'k £ [^l,k ) 



(15) 



Thus, by knowing the dipole transition elements and 
the frequency-components of the electric field of a calibra- 
tion pulse, we can calculate the expansion coefficients of 
the vibrational calibration state created in the electronic 
state a. For convenience, one can merely use a calibra- 
tion pulse that has a constant spectrum over the relevant 
frequencies ufv , giving calibration expansion coefficients 
of the levels according to the values of the dipole tran- 
sition matrix elements, ft is therefore not necessary to 
have detailed information about the calibration pulses, 
but we must know the transition dipole matrix elements 
between the electronic states and a. Fortunately, these 
are typically readily measured, e.g. by fluorescence spec- 
troscopy, or calculated. In contrast, it is more difficult 
to find similar elements between two excited states - ele- 
ments we fortunately do not need to know. 

In the remaining part of this subsection we will be 
considering a single k, and we shall suppress this index 
for clarity. 

For the calibration to be useful, the calibration data 
recorded must enable us to determine the fli's up to a 
common non-zero complex factor. We will accomplish 
this by determining products /3;/3 ; * up to an overall com- 
plex factor, common for all products. The measure- 
ments with pulse 1 gives quantities NiCi(u>)f3i Cp(tu)(3*,, 
precisely allowing us to determine Nipiftp if we know 
Ci(u))C£,(u>). These mapping coefficient products can in 
turn be found, up to a common factor Nc, from the cal- 
ibration data NcQiuj)^ Cp(u)^'* , where the /5f al 's 
are known. By division, we can hence find products of 
the form N'fyftp. To determine the fli's individually, we 
simply set one of these equal to one, use the products 
to determine the rest, and finally normalize. Therefore, 
there is no benefit in knowing one of these products where 
neither I or I' enter in another product. Consequently, we 
must know the range of I's that are populated by pulse 1 
(from its spectral center and width) and be able to find 
products N'fliP*, where all populated Z's enter, and no 
product has both an I and an I' that do not enter in an- 
other. Furthermore, one is free to use data from more 
than one value of u>, as long as each data set can be re- 
lated to data at all other values of u>. This is sufficiently 
fulfilled if one can determine a product NcCi(uj)C*,(uj) 
at one value of u> that contains an I also contained at in 
a product at another value of w. 

In case it is not convenient to populate all these lev- 
els with a single calibration pulse, more can be used in 
subsequent calibrations. Since we can allow for only a 
single undetermined overall complex factor, it is neces- 
sary to relate the data from one calibration pulse to the 
others. For instance, it is sufficient to determine a single 
NciCi(w)C*,(u}) from one calibration CI which can also 
be determined in another calibration C2 to find the ratio 
Nci/N C 2. 



All this having been said, we illustrate the procedure 
with a brief example. Let us say that we are interested 
in reconstructing the vibrational state originating in a 
certain k (which we suppress for clarity), and we desire 
the expansion coefficients (3i with I E {1,2,3,4,5}. We 
employ two calibration pulses CI and C2. CI has a 
constant spectrum in a range sufficient to populate levels 
1-4 and C2 similarly populates levels 3-5. Measuring 
at the two frequencies ui a and u b where other k' 7^ k 
give only negligible contributions, it is found that there 
is good signal strength for the following terms: 

2VciCiMC 2 *(u; a ) ) 
N cl C 2 {u Ja )Cl{u Ja ) } pulse CI 
N cl C 3 {uj b )Cl{uj b ) J 

N C2 C 3 {Lo b )Cl{uj b ) J p 

where we have already divided out the known calibration 
expansion coefficients /3 ; cal /3 ; c , a1 '*, up to a factor absorbed 
in the Nc's- The element Cs(u b )Cl(u) b ) recurring in 
both pulses serves the special purpose of relating the 
two pulses by showing us the ratio of the numbers Nci 
and Nc2- Hereafter, the actual experiment is performed 
with pulse 1. In the table below, we show how one 
obtains the products N'PiPp from measurements and 
calibration data. 



Measurement: 


Used calibration: 


We find: 






N'faft 


N 1 C 2 {uJ a )l32Cl{uJ a )l3l 


iVciC 2 (OC|«) 


N'(3 2 (3* 4 




tv c1 c 3 Mc 4 *m 


N'foPt 




N C2 C 5 (uj b )Cl(uj b ) 


N'(3 3 (3* 4 




N C2 C 3 (Lo b )Ct,(u b ) 





Since we have the freedom of choosing the normaliza- 
tion constant N' , we can arbitrarily set /3x = 1, from 
which we can use the following lines in the table to find 
(in order) N'fa, N'/3 4 , N'{3 S , and from the last line N'/3 5 . 
Finally, the state can be normalized by appropriately 
choosing N' . We have thus found the set of complex 
numbers {kPi}, both their magnitude and phase, and we 
have thereby performed the reconstruction of the state 
formed by pulse 1. 

Having been concerned with how to use small sets of 
measurements, it should be said that in a typical exper- 
imental situation, one will have much more data than is 
minimally required. This abundance of data can be used 
to further improve the statistics of the reconstructed val- 
ues. 

Finally, it may be of note that we can easily find the in- 
dividual mapping coefficients C(w) from the calibration 
data, up to a common factor (rather than just products 
of two coefficients). This is done in the exact same way 
as finding the i.e. by setting one of the Ci(u>) = 1, 

and using the products to find the rest. Thereby dynam- 
ical information about the molecule can be obtained. 
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C. Interference between different k: Case (II) 

We turn our attention to case (II) where not all k can 
be isolated from Eq. (fT2l) by proper choice of ui. This sit- 
uation, encountered in [19] . is illustrated in Fig.[5]for two 
initially populated fc-levels. Here the wave-packet origi- 
nating in k = 1 cannot be found by the above procedure. 

Frequency to 



No pure k = 1 region 



L)q -- 



Both k = and k ■■ 
Interference 



k = dominates 



Figure 5: Schematic illustration of a spectrogram Eq. (|11[) 
in case (II). The vertical lines signify maxima in the signal 
S(w,t)/uj 4 , with contributions from different initial levels ex- 
plicitly shown. In contrast to Fig. 3] we cannot find values of 
w for which only k — 1 contributes. Still, one can first find 
the o/3i's through the procedure described in section II V Bl 
These coefficients can afterwards be used to extract the i/3i 
by using the frequencies in F(f2, to) that arise from optical 
interferences between the k = and k' = 1 signals. If there 
are more overlapping k' signals than the two shown, the yf3i 
from these can be extracted by similar means. 

However, if just a single k can be isolated by choosing 
u> = u>o, we can by the above method find the rj/3z's for this 
k = (usually the ground state) up to a common phase 
factor. In this subsection we describe how the k^ofli can 
be found if the spectrum fulfills certain conditions. For 
transparency, we initially focus on the situation in Fig. [5] 
where we seek the coefficients For our method to 
work, we demand that some of the frequencies f2?',l = 



&k=o k'=v ar i sm g from optical interferences be- 
tween k = and k = 1 , must be spectrally resolved from 
the pure k = and k — 1 signals having f2, ',, — ~' a ' a 



%l> 

Having this requirement fulfilled, we can extract the 
quantities F(U°f,,LUi) = oC;(^i)oA iC;*(u>i) In the 
calibration experiment, we can hence determine the map- 
ping coefficient product qCi (uj\) \C*, (uj\) up to a com- 
mon multiplicative factor and, by using the o/3; found in 
the experiment (using ui — luq), we can find \(3i (using 

Ul = UJl). 

Similarly to case (I), we do not require spectral reso- 
lution of all frequencies and 0^ J, . Indeed, to deter- 
mine all we just need one resolvable for each I 
at some w\, and that the corresponding q/3; is non-zero. 

If more than two values of k overlap, these can be 



found through similar means; either by their interference 
with k = or recursively, i.e. k + 1 is found from k. 
In conclusion, we can completely reconstruct the state in 
case (II). 



D. Case (III) 

The treatments of cases (I) and (II) above were based 
on the assumption that at least a single k can be isolated 
from F(fl,uj) by appropriately choosing ui. If such an ui 
cannot be found, the situation is radically different. Still 
maintaining that we wish to avoid making assumptions 
about the mapping coefficients kCi(uj), one solution is to 
perform several experimental runs at different tempera- 
tures of the initial state. There are two main effects that 
change the spectrograms as the temperature is varied: 
Temperature dependent decoherence processes and ther- 
mal population of the initial state. The most important 
decoherence process will usually be collisional decoher- 
ence (pressure broadening), with a smaller contribution 
from Doppler broadening. In any case, these two effects 
broaden the spectral lines, scaling as VT for fixed par- 
ticle density. For moderate temperature ranges, this is 
only a small effect; indeed we can ignore it altogether as 
long as the required spectral resolution is not washed out. 
In contrast, the initial state populations scale exponen- 
tially with 1/T, greatly changing the spectrogram. It is 
exactly this change we can use to our advantage. Thus, 
the terms in Eq. (fT2|) having k — k' will scale differently 
with temperature for each fc, and differently from terms 
with k 7^ k'. Demanding here that pulse 1 be pertur- 
bative, these variations with temperature simply follows 
the Boltzmann distribution corresponding to the energy 
of the fc'th level in the initial thermal state. By forming 
linear combinations of F(tt, lu) recorded at different tem- 
peratures with weights calculable from the corresponding 
Boltzmann distributions, one can find each of the indi- 
vidual terms with k — k' and with k =/= k' . For example, 
if ki and ki both enters in Eq. (fT2|) . one can find all the 
terms individually; i.e. terms containing (fci, ki), (k2, fo) 
and (fcijfe)- Consequently, it is possible to modify the 
above reconstruction method so it can be used, even if 
all k cannot be selected individually by choosing special 
values of u>. 



V. DISCUSSION 

Having presented above a specific implementation of 
the reconstruction method, we will now discuss how to do 
away with the clarifying, but non-essential, assumptions 
we made in section IIIII We discuss first the effects of 
rotation in subsection IV Al whereafter we will move on 
to discuss the possible generalizations of the method in 
subsection IV Bl 
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A. Rotations 

In the above, we completely neglected rotations, the 
qualitative correctness of which we will now justify. 
Though of general validity, we will for concreteness dis- 
cuss diatomic molecules with rotational energy E rot — 
J( J + 1), where J is the angular momentum quantum 
number. If rotational states up to a certain J ma x are 
appreciably populated, the number of rotational states 
involved is circa {J m ax + l) 2 - In an experiment where 
the temperature is high enough to populate several vibra- 
tional fc states, the number of entering rotational states 
will be large - of the order N ro t ~ 10 4 for the lowest 
vibrational state in the experiment [l9(. All these ro- 
vibrational levels are initially incoherent, and we can 
label each of them after their initial quantum number 
J and its projection M ; The initial state is therefore 
incoherent in all of the indices fc, J and M. Thus, by 
analogy with Eq. ([8]), the dipole moment leading to co- 
herent FW emission only contains contributions of the 

form k,J,Md(t,r) =t(k,J,M'>Pc 3 \T)\d\k,j,M'>Po )t ■ For lin- 
early polarized pulses, each of the one-photon transi- 
tions shown in Fig. [2] changes J by or ±1 while M 
is conserved. Whether AJ = is allowed depends on 
the symmetry of the electronic states — c. In the di- 
atomic case, A J = requires one of the involved states 
in a one-photon transition to have non-zero electronic 
angular momentum. As a consequence, all the tran- 
sition energies between vibrational levels / and I' en- 
tering in Eq. (fTTj) (for clarity, with fc = fc') will be 
split into many & < n a j 1 M \ n, j 2 M \ with Ji and J2 differ- 
ing by at most two. As known from field free align- 
ment, some of these terms will interfere destructively 
after a time 7^ = wh/2BJ and have a full revival af- 
ter T rev — nh/B where J is the mean value of J [25j . 
To elucidate which terms interfere destructively, one can 
write down the frequency dependent dipole moment for 
fixed fc and / quantum numbers in the schematic form 
(i.e. neglecting the difference in populations and map- 
ping coefficients) d ~ Y,j,Ai[ d JMfi + d JM ,+i + djM-i], 
where djM.Aj is the dipole moment contribution orig- 
inating in the state with angular momentum quantum 
number J and having this changed by AJ in the first 
transition. Finding the signal S(ui) oc \d\ 2 /lo a : we find 
N r ot(N rot + 2) ps Nr 0t terms of the form \dj^j\ 2 with 
identical r-dependence; the majority (N^ ot ) coming from 
|dj,Aj=o| 2 - These terms will all stay in phase and give 
rise to a signal at all delay times r, scaling as N^ ot . The 
remaining 8N^ ot — 2N rot ~ 8N^ ot terms will have differ- 
ent r-dependent phases, and will be out of phase a time 
7~di after formation by pulse 1. This makes their contri- 
bution to the signal scale as N rot , except at the fractional 
and full revivals. Summarizing, for delays greater than 
Tdi, but smaller than major fractional revivals, the terms 
with A J = will completely dominate the signal due to 
the scaling N^ ot ^> N rot . Because all these terms have 
the same r-dependence, one can completely ignore the ro- 
tational degrees of freedom. A similar argument can be 



applied to pulse 2 and 3 showing that to avoid diminish- 
ing the signal further due to rotational interference (i.e. 
loose about 8/9 of the signal), it is beneficial to send 
in pulse 2 and 3 so close in time that approximately no 
rotational evolution takes place during their separation 
time. Having accounted for the scalings of the different 
contributions to the signal, we must also account for the 
possible delay times in an experiment. Due to collisional 
decoherence, it is typically only possible to measure the 
spectrogram for values of r much smaller than the rota- 
tional revival time. Conversely, the vibrational periods 
are short enough to go through many cycles before the 
signal is washed away by decoherence. In conclusion, it 
is qualitatively correct to ignore the rotational degrees of 
freedom as long as (i) N rot ^> 1, (ii) A J = transitions 
are allowed, (iii) we limit the r-interval for the Fourier 
transform in Eq. (|12[) so that pulse 1 ends more than 7^ 
before pulse 2 and 3, and (iv) we avoid values of r so that 
pulse 2 arrives at a major fractional rotational revival. 



The rotations also introduces ro-vibrational cou- 
pling, which effectively contributes an additional line- 
broadening. Though usually negligible, this effect scales 
with the absolute temperature as T - more rapidly than 
pressure and doppler broadening. What limits the usable 
temperature is that the required spectral resolution may 
not be washed out by these line broadenings. 



Lastly, it was suggested in section IIVDI to perform 
experiments at different temperatures, altering both the 
vibrational and rotational populations. Being unable to 
resolve the individual rotational states within a vibra- 
tional level, we have taken all these rotational sub-levels 
to be included in their common vibrational level. In this 
way, all mapping coefficients contain sums over rotational 
levels and it is necessary to consider how these mapping 
coefficients change with temperature due to the changed 
rotational content. There could be two causes of such 
a change. Firstly, the transition frequencies & v 'j v i ji in 
Eq. (fTTJjl depend on J. Since A J = 0,±1 in each transi- 
tion, these differences in transition frequencies are rather 
small and may be negligible. If not, one can merely sepa- 
rate pulses 2 and 3 by more than 7^, ensuring that only 
AJ = is measured whereby the transition frequencies 
are independent of J. Secondly, it must be clarified how 
the mapping coefficients and unknown coefficients de- 
pend on the temperature through their rotational states. 
Separating the ro-vibronic states in a product of a rota- 
tional and vibronic part, this is straightforward, but cum- 
bersome. An easy situation arises, however, for linear 
molecules where many rotational states are populated. 
Here, the electronic angular momentum in all states — c 
is of the order 1 and the main contribution to the signal 
comes from J > 1. Then, the transitions depend only 
slowly on J, and the change in difference in rotational dis- 
tributions with temperature may be ignored altogether. 
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B. Generalizations of the method 

Having dealt with rotations, we turn to relaxing the 
other assumptions made in section IHII and the possi- 
ble generalizations of the method. All of these gener- 
alizations are centered around the use of calibration and 
around the linearity in expansion coefficients arising be- 
cause of the perturbative character of the pulses 2 and 3, 
see e.g. Eq. (0). 

First, is easy to see that this structure will be the same 
regardless of whether we are restricted to dipole transi- 
tions or also include higher multi-pole transitions. While 
this may not be of much relevance in four-wave mixing 
experiments, it shows the general principle: While the 
mapping coefficients kCi{uj) become more complicated, 
these are still fully accounted for because of the use of 
calibration pulses. 

Second, the implicit assumption made in Fig. [2] and 
in section IIVI that the photon from pulse 2 enters the 
interaction before the photon from pulse 3 can also be 
relaxed. Similarly to the incorporation of multi-pole 
transitions above, one can sum over the different paths 
through the energy landscape traced out by the differ- 
ent sequences of absorbtion of photons from the pulses 2 
and 3. For instance, if there is temporal overlap of the 
probe pulses 2 and 3, a major contribution to the elec- 
tric field can be caused by an inner product of the form 

(k4' l a\ T )\d\kipi? S> }- Here, |fcVt ) is independent of r, and 
is formed from the initial state by absorption of a photon 
from pulse 2 followed by stimulated emission of a photon 
to pulse 3. Like other energetic paths, this preserves the 
linearity in the kPi's and the signal's dependence on r. 
Indeed, one could even incorporate multi-photon transi- 
tions and still retain the structure of Eq. © and the re- 
construction procedure from section IIVBI provided that 
the truncated perturbation series is an accurate descrip- 
tion at each level of perturbation. If the molecular elec- 
tronic levels are conveniently spaced, one can even use 
another number of pulses for the above procedure, down 
to a minimum of two (one for excitation and one for prob- 
ing) . The special importance of the FW situation is that 
it uses the smallest number of pulses in the widespread 
situation where only two electronic levels = 6 and a = c 
are used. Though our method allows for both different 
numbers of incoming pulses and for multi-photon transi- 
tions, the non-co-linear geometry still plays an important 
role. Its advantage lies not in selecting certain electronic 
states, but rather in avoiding the strong background from 
pulses 1-3. Additionally, it gives a reasonably simple way 
of finding the u-interval that allows us to select a certain 
k' = k in Eq. (fT2j). 

Third, it is unnecessary that the electronic indices 
and a-c denote single electronic states. These could as 
well be groups of states, where the vibrational indices 
k-n simply denotes all states in these groups. With this 
generalization it also becomes apparent that we do not 
need to make the Born-Oppenheimer approximation to 
maintain Eq. @. Without this approximation, we can- 



not speak of vibrational and electronic states, but we 
make instead the more general statement that we find 
the quantum state in the group of energy levels labeled 
by a. 

Fourth, the unknown state does not need to be formed 
by a single laser pulse. What made it convenient to 
use the single laser pulse 1 above, was that the BOX- 
configuration made it possible to distinguish the signal 
from the states having interacted with pulse 1 by the 
direction of the FW-pu\se. Nevertheless, any coherent 
process forming an unknown state in the group of levels 
a will work, provided that there remains population in 
the original group levels (i.e. the states labeled by k), 
and no population is transferred back to these levels by 
the excitation process. 

Fifth, the method above is easily extended to the case 
of correlated initial states, rather than just thermal ones. 
The difference from the above treatment is that there 
can now be quantum correlations between the different 
fc-levels. To treat this situation, one simply writes down 
the density operator corresponding to the vector Eq. (JTJ) , 
retaining the (fc, £/)-correlations. Similarly to the treat- 
ment above, one looks for the coherent FW-teims in the 
emitted radiation with the appropriate direction. If one 
can again find w-values where the mapping coefficients 
are non-zero for a single k only, one can use the method 
above to find the same density matrix elements that we 
found previously. These are the elements arising due to 
the diagonal terms in the energy basis of the initial state. 
The difference, compared to the thermal case, is that 
there is additional information about the state found in 
an cj-region approximately £> fe V from the (k, &)-region. 
This being said, it may well be easier to find the (k, k')- 
correlations in the initial state, instead of trying to ex- 
tract them from the unknown excited state. 

Sixth and finally, there is no requirement that the dif- 
ferent k refers to the same molecule. If one wanted to 
examine a weak emission from a molecule A, it could be 
possible to perform the above experiment on a mixture of 
A with another, more strongly emitting molecule B. In 
the spectrogram region where the signals from A and B 
overlap, the strong electric field from B amplify the sig- 
nal from A rather like the idea of mixing a weak signal 
with that of a strong local oscillator in heterodyne detec- 
tion. The basic principle behind this was demonstrated 
by Engel et al. [2q |. 

By discussing these generalizations, we hope to have 
given the reader the impression that doing quantum state 
reconstruction by perturbative processes with calibration 
is a very general approach with wide applicability. 

In conclusion, we have presented a method to recon- 
struct an unknown vibrational quantum state, coherently 
excited from an initial thermal state. By using one or 
more calibration pulses, we have shown how to circum- 
vent knowledge of the perturbative probing process, and 
we have suggested that this idea may be of much greater 
applicability in quantum state reconstruction. 
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